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In this paper, we present a new experimental progress in brief on the recent observation of the charged 
charmonium-like state Z c (4020) ± states and its iso-spin partner Z c (4020)° in the e + eC -A Tv + ' 0 n~' 0 h c 
process of at the BESIII experiment. The charged Z c (4020) is its decay into 7r ± /i c final state, and carries 
electric charge, thus it contains at least four quarks. The observation of both charge and neutral state makes 
Z c (4020) the first iso-spin triplet Z c state observed in experiment. 

PACS numbers: 14.40.Pq, 13.25.Gv, 13.66.Bc 


I. INTRODUCTION 


In quark model, hadrons are made up of quarks. At the 
lowest configuration, one quark and one anti-quark form me¬ 
son while three quarks do baryon (3]. But hadrons with other 
configurations, called exotic states, are not excluded, such as 
deuteron, which is formed with a proton and a neutron SS- 
The search for these exotic states has stimulated great interest. 
Although no convincible evidence has been found in experi¬ 
mental, strong evidence for mesons that do not fit into the sim¬ 
ple qq-bar scheme of the quark model has been accumulating 
during the past decade, such as X(3872) @-[3], Y(4260) |H|], 
Z( 4430) + HJH i n charmonium region. 

Experimentally, what important to establish exotic states is 
finding clear different signatures compared with the conven¬ 
tional hadrons. One promising method is to search for elec¬ 
trically charged states that can decay into hidden charm final 
states, since such states contain at least four quarks, and there¬ 
fore can not be conventional mesons, if they are mesons. 

In charmonium region, the first such state reported in ex¬ 
periment is Z(4430) + from Belle in 2008 in a study of B —> 
K'K + if{2S) decays @]. This state has not been confirmed 
from BaBar iflltl using the same analysis method. Until the 
result reported by LHCb in 2014 lfl2ll . the Z( 4430) + has been 
confirmed with ambitious significance, and the quantum num¬ 
ber is determined to be 1 + . The first confirmed charged char- 
moniumlike state is Z c { 3900) observed in 7 T^J/ip mass spec¬ 
trum in the e + e _ -A n + n~J/ijj process at BESIII using data 
sample taken at y/s = 4.26 GeV 11131) . This state has been con¬ 
firmed shortly by Belle experiment I14ll and CLEO-c us¬ 
ing same decay mode. Besides, in CLEO-c data, an evidence 
of neutral Z c (3900) has also been reported in the neutral pro¬ 
cess. The observation of Z c (3900) has attracted great interest 
from both theoretical and experimental side. Different mod¬ 
els have been proposed, such as tetraquark consisting of cu 
and cd diquarks fliitl . molecule consisting of D* and D m, 
hadroquarkonium consisting of ud bound to a color-singlet cc 


core 118;], or other configurations [il 9)]. Almost in the same, 
the study of the ( DD *) ± and ( D*D*) ± system in ( DD *) ± 7r 
and ( D*D*) ± n states in the same data sample have been per¬ 
formed, and strong near-threshold peaks has been observed, 
which is named as Z c (3885) ll20ll and Z c { 4025) ll23ll . respec¬ 
tively. 

Apparently, searching for similar states in the invariant 
mass spectrum made up of one pion and other charmonium 
states such as h c , XcJ, and so on is interesting. Any ob¬ 
servation may provide more information on understanding 
these states, and helpful in building relationship between those 
states. Such studies have been performed at BESIII in both 
e + e _ -A TT + n~h c and n°TT°h c processes l2ll[22ll . A distinc¬ 
tive structure, named as Z c (4020) has been observed in the in¬ 
variant mass spectrum of 7 rh c . The mass of Z c { 4020) is close 
to that of the mentioned Z c { 4025) ± observed in (D*D*) ± tt 
states. Here we just introduce the observation of both electri¬ 
cally charged and neutral Z c ( 4020) at BESIII. 
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FIG. 1. The comparison between <j(e + e -A 7r + 7r h c ) and 
a{e + e~ -A 7r + 7r ~ J/ip). The errors are statistical only. 














2 


TABLE I. The cross section of e + e“ —> mvh c at different energy points and the ratios (R nnhc ) between charged and neutral modes. For the 
first three energy points, besides the upper limits, the central values and the statistical errors are also listed. The second errors are systematic 
errors and the third ones are from the measurement of branching fraction of h c -A 7 r] c {2 61 . 


y/S( GeV) <r(e + e -A 7r + 7r h c ) 

<r(e 

+e 

r -7 

7r°7T 
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c) 

-^7T7r/l ( 



3.900 

0.0 ±6.0 or <8.3 



- 




- 



4.009 

1.9 ± 1.9 or <5.0 



- 




- 



4.090 

0.0 ±7.4 or <13 



- 




- 



4.190 

17.7 ±9.8 ±1.6 ±2.8 



- 




- 



4.190 

17.7 ±9.8 ±1.6 ±2.8 



- 




- 



4.210 

34.8 ±9.5 ±3.2 ±5.5 



- 




- 



4.220 

41.9 ± 10.7 ±3.8 ±6.6 



- 




- 



4.230 

50.2 ±2.7 ±4.6 ±7.9 

25.6 

± 

4.8 ± 

2.6 

± 

4.0 0. 

54 ±0.11 

± 

0.06 

4.245 

32.7 ± 10.3 ±3.0 ±5.1 



- 




- 



4.260 

41.0 ±2.8 ±3.7 ±6.4 

24.4 

± 

5.2 ± 

3.2 

± 

3.8 0. 

63 ±0.14 

± 

0.10 

4.310 

41.0 ±2.8 ±3.7 ±6.4 



- 




- 



4.360 

52.8 ±3.7 ±4.8 ±8.2 

36.2 

± 

6.5 ± 

4.1 

± 

5.7 0. 

73 ±0.14 

± 

0.10 

4.390 

41.8 ±10.8 ±3.8 ±6.6 



- 




- 



4.420 

49.4 ± 12.4 ±4.5 ± 7.6 



- 




- 




II. MEASUREMENT OF e+e" -A n + n~h c AND 

OBSERVATION OF e + e“ -A n°n°h c 

The process of e + e _ -A n + iT~h c was first observed by 
CLEO-c experiment using 586 pb ~ 1 data sample at center-of- 
mass (CM) energy y/s = 4.170 GeV. Together with the scan 
sample around 4.000 GeV and 4.260 GeV, the cross section 
shows a hint of rising at 4.260 GeV tf25ll . as shown in Fig¬ 
ure [I] An improved measurement may shed light on under¬ 
standing the nature of the Y states and also used to exam pos¬ 
sible charged intermediate charmoniumlike states. 

The BESIII experiment, works at the Beijing Electron 
Positron Collider (BEPC) II, is a tau-charm factory which can 
collect data from yfs = 2 GeV to 4.6 GeV (2M . Starting 
from the end of 2012 to 2013, BESIII detector collected about 
3.3 fb _1 data samples at 13 CM energies from 3.900 to 4.420 
GeV, which are listed in Table I. 

Based on these data samples, the process of e + e _ —> 
nnh c has been studied using full reconstruction of the 
final state particles. Here the h c is reconstructed via 
its electric-dipole (El) transition h c —> 7 ?y c with rj c 
reconstructed with its 16 hadronic decay modes: pp, 
TT + TT~K + K~, TT + 'K~pp, 2 (K + K~), 2(7T + 7T - ), 3(7T + 7T _ ), 
2(tT + TT~)K + K~, KsK^TT^, KsK ± TT ZfZ TT + TT~, K + K~TT°, 
K + K~p, pp 7 T°, n + TT~ri, 7 r + 7 r _ 7 r° 7 r°, 2('K + ir~)i r i, and 
2 ( 7 r + 7 r _ 7 T°). Here Kg is reconstructed in its 7 r + 7 r~ decays, 
and t ] in its 77 final state. 

Figure[2}a) shows as an example the scatter plot of the mass 
of the r] c candidate versus that of the h c candidate at the CM 
energy of 4.260 GeV (left panel), as well as the projection 
of the invariant mass distribution of 777 in the rj c signal re¬ 
gion (right panel) in e + e~ —>- TT + n~h c process, where a clear 


signal is observed & The 7 r] c mass spectrum is used to ex¬ 
tract the number of Ti + TT~h c signal events by performing a fit. 
The fit to the 4.260 GeV data is also shown in Figure [2] At 
the energy points with large statistics (4.230,4.260, and 4.360 
GeV), the fit is applied to the 16 r] c decay modes simultane¬ 
ously, while at other energy points, the fit is performed to the 
mass spectrum summed over all the r] c decay modes. With 
the fitted number of hc events and the corresponding detect 
efficiency as well as the integrated luminosity at each energy 
point, the cross section of e + e _ -A n + TT~h c at different en¬ 
ergy points are listed in Table [Q 

The cross section of e + e _ —> Tr°TT°h c is also studied using 
3 data samples at CM energy of 4.23, 4.26 and 4.36 GeV. Fig¬ 
ure [ 2 a) shows the scatter plot of the rj c candidate versus that 
of the h c candidate summed over 3 energy points, and Fig- 
ure[ 2 b) shows the projection of the invariant mass distribution 
of 'yrjc in the rj c signal region, clear h c signal is also observed 
in the neutral process. The 777 ,. mass spectrum summed over 
all the T] c decay modes is fitted simultaneously at the 3 energy 
points to extract the number of ir°ir°h c signal events. The 
cross sections of e + e _ TT°n°h c are also listed in TableQ] 

The comparison of the cross section of e + e _ —> Tr + TT~h c 
measured at BESIII with the cross section of e + e —> 
TT + n~J/ , ip measured at Belle experiment using initial state 
radiative (ISR) method fl4ll is shown in Figurt|T| The cross 
section of e + e _ —> n 0 TT 0 h c is at the same order of magni¬ 
tude of 7 r + 7 r _ J/ip but with a different line shape. There is a 
possible broad structure at high energy with a local maximum 
around 4.23 GeV. Since the Y(4260) was established from the 
cross section of e + e _ -A J/ijj, the different line shape 

e + e _ -A 7r + 7r ~h c makes the understanding of Y states more 
complicate. 
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(a) 


(b) 


FIG. 2. The scatter plot of the mass of the r; c candidate versus that of the h c candidate at the CM energy of 4.260 GeV (a) and The projection 
of the invariant mass distribution of 7 r] c in the rj c signal region, the dots with error bars are data, and the solid curve shows the best fit result 
(b). 
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FIG. 3. The scatter plot of the mass of the ?; c candidate versus that of the h c candidate (a) and the M™o°o distribution for the events with an 
rjc candidate(b). 


Combining the Bom cross sections of e + e” —> TT 0 n°h c 
with that of charged mode, the ratios of cross section between 
charged and neutral mode () at each energy points are 
calculated, as listed in TableQ] In the calculation, the common 
systematic uncertainties in the two measurements have been 
canceled. The combined ratio for the 3 energy points , 

obtained with a weighted least squares method lt27ll . is deter¬ 
mined to be 0.63 ± 0.09, which agrees with the expectation 
based on isospin symmetry within a systematic uncertainty of 
0.5. 


III. OBSERVATION OF THE CHARMONIUMLIKE STATE 

Z c (4020) 

At the three energy points with large statistics, the in¬ 
termediate states are studied by examining the Dalitz plot 


of the selected hc candidate events. Figure |4ja) shows the 
projection of the M(7 t ± /i c ) (two entries per event) distribu¬ 
tion for the T: + 'K~h c candidate events summed over 3 energy 
points. There is a significant peak at around 4.02 GeV/c 2 
(named Z c (4020) ± hereafter), and there are also some events 
around 3.9 GeV/c 2 (inset of Figure S]a)), which may be the 
candidate events of the observed Z^^OO) 7 * 1 in e + e - —> 
tt + 7t~ J/ip a . In the 7T7T mass spectrum, there is no obvious 
structure. 

Similarly, in the M£® coll | max (one entry per event, the larger 
one in two % c candidates is retained) mass distribution for 
the e + e _ —> ir 0 ir°h c events, there is also a significant peak, 
which could be the isospin partner of Z c (4020) , at around 
4.02 GeV/c 2 (named 2' c (4020) 0 hereafter). 

M(ir^hc) or AFo coll |max distributions summed over the 16 
rj c decay modes is fitted simultaneously at the three energy 
points to extract the parameters of the structures and the sig- 
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FIG. 4. Sum of the simultaneous fit to the M 7r ± hc distribution (a) and sum of the simultaneous fit to the M^ coll |max distribution at yfs = 
4.23, 4.26 and 4.36 GeV as described in the text(b). 


TABLE II. Cross sections of e + e — > 7tZ c (4020) — > mvh c measured from the BESIII experiment and ratios R^Zcf^o) mentioned in the 
text. The first errors are statistical errors, the second errors are systematic. Here, cr® TZ ( 4020 )* denotes a B (e + e~ — > Z c (4020) ± —> 
'K + n~h c ) and °’^o Zc (4020) 0 denotes a B (e + e~ —> n° Z c (4020)° —> 7r°7r°h c ). 


v/i(GeV) 

obs 

,l Z c (4020)± 

Z c (4020)^ (pIO 

obs 

n Z c ( 4020)° 

Z c (4020) 0 

f?7rZ c (4020) 

4.23 

114 ±25 

8.4 ±1.8 ±2.6 

22 ±7 

6.5 ±2.2 ±0.7 

0.77 ±0.31 ±0.25 

4.26 

76 ± 17 

7.0 ±1.6 ±2.0 

23 ±8 

8.5 ±2.9 ±1.1 

1.21 ±0.50 ±0.38 

4.36 

67 ± 15 

9.8 ±2.2 ±2.9 

17 ±7 

9.9 ±4.1 ±1.3 

1.00 ±0.48 ±0.32 


nal events yield.The mass and width of Z c (4020) ± are mea¬ 
sured to be 4023.6 ± 2.2 ± 3.8 MeV/c 2 , and 7.9 ± 2.7 ± 2.6 
MeV, respectively, with statistical significance of 8.9 (see Fig¬ 
ure [4j a)), while the mass of Z c (4020)° is determined to be 
(4023.6 ± 2.2 ± 3.8) MeV/c 2 with a statistical significance 
larger than 5cr (see Figure®!))). Here, the width of Z c (4020)° 
is fixed as its charged isospin partner due to the limitation of 
signal events. The detail information about the systematic un¬ 
certainties in mass and width measurement of Z c (4020) can 
be found in published paper lElial]. In the charged mode, 
the possible decay of (3900)^ —> n^hc has been explored 
by adding a Z c (3900) structure in the fit with mass and width 
fixed to the BESIII measurement @|. It is found the statisti¬ 
cal significance of Z c ( 3900) is less than 2.1<r (see the inset of 
Figure 0] a)) 0. 


The signal yields of Z c (4020), as well as the measured 
cross sections for e + e _ — > irZ c (4020) —> Tnrh c and the ra¬ 
tios IZttZc between neutral and charged modes at three en¬ 
ergy points, are listed in Table [IT] The common systematic 
uncertainty in the ratio calculation has been canceled. The 
combined ratio TZ n z c for 3 energy points is determined to 
be (0.99 ± 0.31) with the same method as for the com¬ 
bined IZ^TThc’ which is consistent well with the expectation 
of isospin symmetry = 0.5 within 1.6er. 


TABLE III. The resonance parameters of the Z c (4020) and 
Z c (4025). 

State Mass (MeV/c 2 ) Width (MeV) 
.Z' c (4020) ± 4022.9 ± 0.8 ± 2.7 7.9 ± 2.7 ± 2.6 
Z c (4020)° 4023.9 ± 2.2 ± 3.8 Fixed (=7.9 MeV) 

Z c { 4025) ± 4026.3 ±2.6 ±3.7 24.6 ± 5.6 ± 3.7 


IV. SUMMARY AND OUTLOOK 

In summary, the charmoniumlike structure, Z c { 4020), has 
been observed in irh c system for both charged and neutral 
modes. It is not only a multi-quark state, but also a isospin 
triplet. Table [III] shows the measured resonance parameters 
of the charged and neutral 2 C (4020) state. The mass dif¬ 
ference between neutral and charged 2).(4020) is 1.0 ± 2.2 
(stat.) MeV/c 2 , which in good agreement with zero within un¬ 
certainty. Furthermore, the measured Born cross sections for 
<r(e + e _ -A n°n°h c ) are about half of those for a(e + e~ —>■ 
TT + TT~h c ) within errors, while the measured Born cross sec¬ 
tions <r(e + e“ —> 7r°7r°/i c ) for charged and neutral modes. 

Furthermore, the measured Born cross sections for 
cr(e + e~ — > n 0 n°h c ) are about half of those for a(e + e~ -A 
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Tr + TT C h c ) within errors, while the measured Born cross sec¬ 
tions <r(e + e - —i 7 tZ c (4020)) for charged and neutral modes 
are generally consistent. It indicates that there is no obvious 
isospin violations in nnh c and 7rZ c (4020) systems. 

The observation of charged Z c ( 4020) confirms that there 
is an extra multi-quark state except the observed Z c (3900). 
In particular, c is the first observed neutral charmoniumlike Z 
state, and it is helpful to understanding the nature of Z c and 
even other XYZ particles. For instance, the observation of 
theoretical expectations of Z c (4020) and X(3872) combined 
with Y (4260) transition to X(3872), which establish contact 
with XYZ particles, issue the challenge to interpret these par¬ 
ticles theoretically. 

Compared to the ZJ 4025) observed in the process 
e + e~ —> (D*D*) ± tt^ : [l23j], the resonance parameters of 
Z c (4020) lEMi agree with Z c (4025) (H with 1.5cr. If 
they are the same state, the coupling of Z c (4020) to D*D* 
final state is about 12 times larger than to nh c final state. To 
clarify whether they are the same particle, one needs a further 


sophisticated analysis with a coupled channel technique. 

The observation of charged charmonium-like family mem¬ 
bers together with the similar observation in bottomonium 
family provides important evidence for the existence of 
hadrons beyond the conventional quark model and QCD the¬ 
ory, and helps in the improvement of our understanding about 
the fundamental structure of the matter. 

In the near future, more data between 4.0 and 4.6 GeV at 
BESIII experiment |24| will be accumulated for further study 
of XYZ states. The property of Z c (4020) will be studied 
more carefully. Meantime, the line-shape of cross section of 
e + e~ — i TTirh c and e + e~ —> 7rZ c (4020) will be investigated, 
which can shed light on the understanding of this state and 
even other XYZ particles. 
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